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Abstract 

We present a study of the LHC sensitivity to the W boson mass based on simulation 
studies. We find that both experimental and phenomenological sources of systematic 
uncertainties can be strongly constrained with Z measurements: the lineshape, doz/dm, 
is robustly predicted, and its analysis provides an accurate measurement of the detector 
resolution and absolute scale, while the differential cross-section analysis, d^Oz/dydpT, 
absorbs the strong interaction uncertainties. A sensitivity Sm^ 7 MeV for each decay 
channel {W ev, W —>■ jJ-V), and for an integrated luminosity of 10 fb^', appeai^s as a 
reasonable goal. 



1 Introduction 



The Standard Model (SM), now computed at two-loop precision [1,2], is a very predictive framework. 
Its most precisely measured parameters aQED, G^, and mz, provide constraints on the W boson and 
top quark masses, which can be confronted with measurement. Injecting the measured value of the W 
mass and the measured Z boson couplings, a definite prediction is given for the top quark mass [3]. 
This prediction, together with the discovery of the top quai^k at a compatible mass [4,5], has been a 
major achievement in high energy physics. 

The measured values of the W boson and top quark masses are now more precise than their quantum 
predictions, and provide non-trivial constraints on the gauge symmetry breaking sector. In the SM, 
this translates into limits on the Higgs boson mass [6]. Beyond the SM, constraints are given on the 
contributions of other heavy particles, like supersymmetric particles [7]. 

The W mass has been measured at UA2 [8], LEP [6], and the Tevatron [9]. The recent measure- 
ment by the CDF Collaboration gives mw = 80.413 ± 0.048 GeV, yielding a current world average of 
mw = 80.398 ±0.025 GeV [10]. In the SM, the resulting Higgs boson mass uncertainty is about 50%. 
Any further improvement in this measurement will translate into more precise indirect predictions. 

The present paper discusses the LHC prospects for the W mass measurement. The expected W 
cross-section at the LHC is about 20 nb [11]. In 10 fb^^ of data, a benchmark for one year of in- 
tegrated luminosity during the first years of stable running, around 4 x 10^ W events will be selected 
in each exploitable decay channel (W — > ev,lJ.v), providing a combined statistical sensitivity of about 

1 MeV. Previous estimates [12-14] of the systematic uncertainties affecting this measurement how- 
ever amount to 5mw '-^ 20 MeV per experiment, and to a combined uncertainty of 5mw ~ 15 MeV. 
The main sources are the imperfect determination of the experiments absolute energy scale, and the 
uncertainties in the W boson kinematical distributions (rapidity, transverse momentum), which in turn 
stem from proton structure function uncertainties and higher orders QCD effects. 

The purpose of this paper is to re-investigate the possibilities to measure the W mass with the great- 
est possible precision. As is known from the Tevatron, the uncertainties can be significantly reduced 
using Z boson measurements; this approach will be employed here, with modifications and improve- 
ments suggested by the high Z statistics expected. Although our discussion is general, most of our 
arguments rely on the expected performance of the ATLAS experiment [15]. 

The paper is structured as follows. Section |2] summarizes the W mass fitting procedure, lists the in- 
gredients needed to describe the W distributions used in the fit, and gives a general description of how 
these ingredients can be determined. The sources of uncertainty are then discussed in turn, in Sec- 
tion [3] (experimental uncertainties), Section|4](theoretical uncertainties), and Section |5] (backgrounds, 
underlying event, and effects related to the machine operation). Correlations between these effects are 
discussed in Section [6l and the results are given in Section |7] Section [8] concludes the paper. 

2 General discussion 

This section discusses our technical set-up, the W and Z event selection, the mass fitting procedure, 
and the problem of controlling all ingredients entering in the definition of the fitted distributions. 
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2.1 W and Z production. Event generation and simulation 



Throughout this paper, W and Z boson samples, and their distributions and acceptances are computed 
using the PYTHIA general purpose event generator [16]. On top of PYTHIA, the treatment of pho- 
ton radiation in W and Z decays is done via an interface to PHOTOS [17]. The size of the expected 
samples are computed assuming the NLO W and Z cross-sections, as obtained from RESBOS [18]. 
These choices are not unique, and the simulation of physics processes at the LHC, in particular non- 
perturbative strong interaction parameters, will obviously need to be adjusted using the forthcoming 
data. In this analysis, the effects of the corresponding uncertainties are estimated either by changing 
parameters in these programs, or by distorting the output distributions according to our assumptions. 

When referring to "fast simulation", we mean a simplified simulation of the ATLAS detector response 
using scale factors and Gaussian resolution functions, applied to the generator-level information ob- 
tained above [19]. When referring to "full simulation", we mean the complete simulation of the 
ATLAS detector using GEANT4 [20]. In our discussions below, and in the absence of real physics data, 
we often treat our fully simulated event samples as data samples, and the fast simulation samples as 
their Monte-Carlo simulation. The different detector response in fast and full simulation allows to 
emulate the realistic situation where the imperfect detector simulation is adjusted during data taking. 

2.2 Signal selection and fitting procedure 

At hadron colliders, W and Z events can be detected and reconstructed in the eVg, /iv^, ee, and HjU. 
final states. The hadronic modes suffer prohibitively large background from jet production; T modes 
can be detected but the T-lepton decay produces additional undetected particles in the final state, di- 
luting the information that can be extracted from these modes. In W events, the observables most 
sensitive to mw are: 

• The reconstructed lepton transverse momentum, pj ; 



The transverse momentum of the neutrino, p^, is inferred from the transverse energy imbalance, cal- 
culated from a summation of energy in all calorimeter cells. Electrons are measured using the inner 
detector (ID) and electromagnetic calorimeter (EMC). They are reconstructed and identified with an 
efficiency of about 65%, while rejecting background from jets up to one part in 10^; in W decays, the 
energy resolution is about 1.5%. For muons, the ID is used together with the muon spectrometer; the 
reconstruction efficiency is about 95% and the relative momentum resolution about 2% [21]. 

The W signal is extracted by selecting events with one reconstructed isolated, high-pj- lepton (electron 
or muon), large missing transverse energy (due to the undetected neutrino), and low hadronic activity. 
In the following, we require p^>20 GeV, |T]f | < 2.5, £'™'">20 GeV, and require the hadronic recoil 
(defined as the vector sum of all calorimetric transverse energy opposite to the reconstructed W decay 
products) to be smaller than 30 GeV. These selections have a total efficiency (trigger and selection) of 
about 20%, providing a sample of about 4 x 10^ events in each decay channel. The backgrounds are 
at the percent level. Table [T] summarizes these numbers. The pj and ni^ distributions obtained with 
fast simulation after the W event selection are shown in Figure [T] 

Based on these distributions, niw can be extracted by comparing the data to a set of models (or template 
distributions) obtained from W event generation followed by a fast simulation of the decay particles. 



• The reconstructed W transverse mass, ntj = 
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Channel W ^ Iv Z ^ II 

Cross-section (pb) 19800 1870 

Lepton T] acceptance 0.63 0.51 

Selection eff. ~ 0.2 ~ 0.2 
(including acceptance) 

Expected statistics (10 fb"^) 4 x lO' 3.5 x 10^ 

Table 1: Cross-section, T] acceptance, total selection efficiency (averaged for electrons and muons) 
and expected sample size for 10 fb^\ in each decay channel. 




Figure 1: Distributions of the lepton transverse momentum, pj, and of the W transverse mass, nij, 
after typical W event selections (cf. text). The Jacobian edges in these distributions provide sensitivity 
to the W mass. 
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Figure 2: value as function of the tested value of mw- Each dot represents a comparison between 
the data and the template distribution obtained for a given mw. The curve is the fitted parabola. 



The different template distributions are obtained by varying the value of the W boson mass parameter 
in the event generation. The statistical comparison of the data to the templates can be performed in 
various ways; throughout this study we will use a simple binned test. The x"^ quantifying the 
compatibility of a given template distribution with the data is defined as follows: 

( ob, _ exp.2 
i i 

where n^^ and n°^^ are the number of expected and observed events (in the template distribution and 
in the data, respectively) in bin / of the pj or spectrum, a,- is the expected resolution, and the sum 
extends over all bins in the fitting window. The Gaussian approximation used above is justified for 
large statistics, which is the case we consider here. 

After all x^ evaluations, a parabola is fitted through the x^ values as a function of m^/. The procedure 
is illustrated in Figure |2l With the statistics given in Table [T] each channel provides a statistical preci- 
sion of about 2 MeV for data corresponding to an integrated luminosity of 10 fb^^. 



2.3 Required inputs 

For the above procedure to work in practice, one must predict the pj and distributions as a func- 
tion of the W mass. These distributions are however affected by many effects, which need to be 
included correctly in order to avoid biases in the mass fit. The needed inputs are listed below. 

• Experimental inputs: the energy scale and resolution need to be known in order to describe the 
Jacobian peak coiTcctly (position and spread). Electron and muon reconstruction efficiency effects 
also distort the spectra, if this efficiency is pj and T] dependent. 



• Theoretical inputs: the W rapidity distribution, y , affects the and pj distributions. The trans- 
verse momentum of the W, pj, directly affects the pj spectrum; its impact is weaker on the 
spectrum. The and p^ distributions depend on the proton structure functions and on higher-order 
QCD effects. The lepton angular- distribution in the W rest frame is of importance for both pj and 
my and changes with the W polarization [22]. Finally, QED effects (photon radiation in the W decay) 
shifts the lepton pj downwards. Since the radiated photons are mostly collinear to the charged decay 
lepton, the impact on electrons and muons is different: the measured muon momentum entirely re- 
flects the momentum loss by radiation, whereas the electron energy, measured essentially in the EMC, 
includes most of the radiated energy. 

• Environmental inputs: these include, among others, backgrounds surviving the W selection, un- 
derlying event and pile-up effects on reconstructed energies and momenta, random neutron hits in the 
muon spectrometer ("cavern background"), and the impact of a non-zero beam crossing angle. In all 
cases, imperfect modelUng of these inputs biases the event reconstruction, leading to distorted pj and 

distributions. 



2.4 Propagation of systematic uncertainties 

The impact of underlying physics mechanisms affecting the W mass determination is estimated by 
producing template distributions of pj and rn^ unaware of the effect under consideration, and fitting 
them to pseudo-data including this effect. The resulting bias (i.e. the difference between the injected 
and fitted values of miy) gives the corresponding systematic uncertainty. 

In the simplest case, a physics effect affecting the distributions (for a given value of the W mass) can 
be summarized by a single parameter. In this case, the induced systematic uncertainty is simply given 
by: 

dmw = ^^{dreia) (2) 

OrelOC 

where a is the parameter controlling the parasitic physics effect, 5reiOC its relative uncertainty, and 
5mw the induced systematic uncertainty on the W mass. When applicable, we will quote the un- 
certainty SreiCC, the derivative dmw/dreiCC and the estimated 5mw- As a convention, we normalize 
dmw/dreiCC in MeV/%. 

Sometimes, however, a single parameter is not sufficient. The uncertainty 5mw is then the result of 
all parameter uncertainties and their congelations: 

Omw =2^T ^ {OrelCCi){Oreiaj)pij. (3) 

This happens when the systematic is parametrized by a (sometimes empirical) function. In this case, 
we choose to quantify the impact by Monte-Carlo propagation: we generate random configurations of 
the ttj, within their uncertainties, and preserving their correlations; for each configuration, we produce 
the corresponding pseudo-data, and fit them to the unaffected templates. The spread of the distribution 
of the fitted mw values gives the contribution to 5mw- 



6 



2.5 The impact of Z boson measurements 

The LHC will produce a large number of Z events. Their selection is rather straightforward, requir- 
ing two reconstructed isolated, high-pt leptons {pj>20 GeV, \ric\ < 2.5), and low hadronic activity 
(hadronic recoil smaller than 30 GeV). 

For each useful decay mode (Z — > ee,/i/i) and for ~ 10 fb^\ around 3.5 x 10^ events should survive 
selections. This represents a factor 10 less than the expected W statistics, but the fact that Z events 
are fully reconstructed largely compensates this deficit. Cross-sections and statistics are summarized 
in Tabled 

The precise knowledge of the Z mass and width will allow to determine the lepton energy scale and 
resolution precisely. Exploiting the energy distribution from the decay leptons will also allow to de- 
termine the scale's energy dependence (i.e, the linearity of the detector response), and the energy 
dependent resolution function. Once this is achieved, the Z transverse momentum will also serve to 
scale the measured hadronic recoil to the Z; together with the measured lepton transverse momentum, 
this defines the missing transverse energy. Finally, "tag and probe" methods [23] will allow to deter- 
mine the lepton reconstruction efficiency. 

Although most of the QCD mechanisms affecting W distributions carry significant uncertainty [24], 
they affect W and Z events in a similar way. This is the case for non-perturbative contributions to 
the W transverse momentum distributions, but also for parton density (PDF) effects: at the LHC, the 
W and the Z are essentially sensitive to high-2^ sea partons, and a variation of these parameters will 
affect the W and Z distributions (in particular , y^) in a correlated way. Hence, the measurement of 
the Z distributions will help to control the W ones. 

The simulation of QED radiation in W and Z decays was much improved recently [17, 25]. Still, 
the measurement of this process (through e.g. Z — > tij) will allow to confirm the predictions. Other 
sources of uncertainty (e.g. backgrounds and underlying event) will also be controlled by auxiliary 
measurements at the LHC. 

The following sections attempt to quantify the above arguments. 

3 Experimental uncertainties 

This section assesses the effect of efficiency and resolution in the reconstruction of leptons and missing 
transverse energy. 

3.1 Lepton scale and resolution 

The Z boson resonance has been measured very precisely at the lepton colliders during the 90's [3]. 
The Z boson mass and width can be exploited as an absolute reference to determine as precisely as 
possible the detector energy scale, its linearity and resolution. 

The basic method is rather simple, and consists in comparing the position and width of the observed 
mass peak in reconstructed dilepton events with the Z boson parameters. A shift of the observed po- 
sition of the mass peak, with respect to the nominal Z peak position, is corrected for by scaling the 
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detector response, hence determining the detector absolute scale; the additional spread of the mass 
distribution, as compared to the natural Z boson width, is used to estimate the resolution. 

The high statistics expected at the LHC, however, imposes a number of refinements. First, the scale 
obtained as above is averaged over the lepton kinematical spectrum, whereas an energy-dependent 
scale is needed for a correct description of the Jacobian distributions in W events. Secondly, lepton 
energy resolution effects induce a small but non-negligible shift in the di-lepton invariant mass distri- 
bution. This shift needs to be subtracted before converting the scale measured from the Z invariant 
mass distribution into the scale used to describe the Jacobian distributions in W events. The resulting 
method has been described in detail in Ref. [26], and is summarized below. 

3.1.1 Average detector scale 

We first illustrate the energy-independent method, providing an average detector scale. Using the 
PYTHIA event generator [16], we produce a set of template histograms corresponding to generator- 
level Z lineshapes. The decay leptons are smeared and decalibrated with different energy scale factors 
a and resolution functions a. For definiteness, we consider calorimeter-like resolution functions 
parametrized as o{E) = a x \fE. At this stage, a is independent of the lepton energy. These tem- 
plates are to be compared to the data; for our tests, we use an independently simulated sample as 
pseudo-data. 

A test is then performed between the pseudo-data and each of the template histograms, as in Sec- 
tion |T2l This results in a two-dimensional grid as a function of the smearing parameters. At the 
vicinity of the minimum, a paraboloid can be fitted through the points, and the parameters of this 
pai^aboloid give the estimates of the true values of a and a. 

The method is tested on a fully simulated Z — > sample, corresponding to 30700 events with 
85 < itiee < 97 GeV; the mass resolution can be treated as Gaussian over this range. We find an 
average resolution parameter a = 0.142 ±0.003, and an average mass scale a = 1.0038 ±0.0002. 
Figure |3] illustrates the result, where the fully simulated Z ^ ee lineshape is compared to an example 
template histogram assuming a = I and a = 0. 12, and to the best fit result. Very good agreement is 
obtained; moreover, the "measured" scale and resolution parameters coincide with the values found 
when comparing the reconstructed electron energies to their generation-level values; the electron cal- 
ibration in the fully simulated sample underestimates the true energy by 0.4%. 

Assuming an inclusive Z production cross-section of 2 nb per leptonic decay channel and an integrated 
luminosity of 10 fb^^ the average scale and resolution parameters can be controlled with a relative 
precision of SreiCC = 2 x 10^^ and dreia = 2 x 10^^. Note that these values are not far from the actual 
uncertainty of the Z boson parameters. As far as the absolute scale is concerned, a coiTclation between 
the induced W-mass systematic uncertainty and the Z boson mass uncertainty might finally appear. 

As discussed in the introduction to this section, the method illustrated here has an important shortcom- 
ing: it only provides a scale averaged over the pj- distribution expected in Z events, which differs from 
that expected in W events. The averaged scale is applicable to W events only in the absence of any 
non-linearity in the detector response. In order to correctly propagate the Z calibration measurement 
to the W sample, the scale thus needs to be measured as a function of energy. This is discussed next. 
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3.1.2 Linearity: energy dependent scale and resolution 

The above method can be extended as follows. The data and the templates are classified as a function 
of the lepton energies. This leads to templates and pseudo-data labeled corresponding to the 
event categories where one lepton falls in bin /, and the other in bin j. The scale factor a,y and the 
resolution parameter a,y are then fitted in every bin. 

In case of small non-lineaiities of the calorimetry response (i.e. a,-y,a,,ay very close to 1), we can 
then derive the a, from the a,y, writing in first order approximation that the mass peak decalibration 
results from the decay lepton decalibrations: 



a,'ymi2 = y 2 UiE I a jE2 {I — COS 612), or (4) 
aij ~ (a; + ay)/2 (5) 

Writing this for every gives a linear system which can be solved using least squai^es. 

As for the resolution, the following linear system holds, neglecting the small contribution from the 
angular terms in the expression of the invariant mass resolution: 

^ ^ (6) 

which can again be solved using least squai^es, yielding the a,. We thus obtain the energy-dependent 
resolution function, independently of the form used to produce the templates. 
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Figure 4: Left: reconstructed absolute energy scale, as a function of energy (points with error bars). 
The fuU line gives the injected function, representing the effect of non simulated passive material in 
front of the calorimeter. The dashed line is an empirical function fitted through the points; the dot- 
dashed line shows the result of an energy-independent analysis, missing the non-linearities. Right: 
reconstructed energy resolution (points with eiTor bars). The full line is the true resolution function, 
of the form a{E)/E = a/\/E + b; the dashed line is the reconstructed function. The dot-dashed Une, 
assuming no constant term {b = 0), is strongly excluded. 



Examples of results that can be achieved are shown on Figure |4l With energy bins defined as intervals 
of 5 GeV, and a integrated luminosity of 10 fb^^ the scale parameters are reconstructed with a preci- 
sion of 2 X 10^"*, as estimated from the RMS of the a residuals with respect to the injected function. 
Similarly, the resolution parameters are reconstructed with a precision of 2 x 10^^. 

3.1.3 Propagation to my/- Smwiccc), 5mw{oe) 

Assuming that bin-to-bin variations of the scale do occur with a spread of 2 x 10^^, we can compute 
the impact of such variations on the measurement of mw- 

As described in Section |2l we perform a set of toy measurements, using the electron transverse mo- 
mentum as observable, templates with varying mw values but with a perfectly linear scale, and pseudo- 
data with fixed mw, but containing non-linearities. 

First of all, we can study the mw bias as a function of the error on the average absolute scale. Not 
surprisingly, we find a strong dependence: 

dmw 



drelCCC 

as illustrated in Figure [5] 



800 MeV/%, 



In the case of an energy-dependent scale, the uncertainty on mw is obtained by injecting random, 
energy-dependent decalibrations in the pseudo-data, with a spread corresponding to the result of the 
analysis of Section [3.1.2l With 480 independent exercises of this type, we obtain a distribution of m^' 
as shown on Figure [6l The scale-induced W mass uncertainty is given by the spread of this distribu- 
tion, and is 5mw{oCi) = 4 MeV. 
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8a, (%) 5C| (%) 

Figure 5: Left: bias on mw, — m^"'", as a function of the relative bias on ai, da^ = [af" — 
a'/"") / af"'' . Right: bias on as function of the resolution bias, 5oi = (a/" — a'/"'') / a'/"'' . A 
Unear dependence is observed in each case, with dmw / dreiCCe = 800MeV/% and dmw/dreiOi = 
0.8 MeV/%. 




Figure 6: Distribution of w?^ , for 480 exercises with energy-dependent scale parameters randomly 
drawn within theii" uncertainties. The true mass is mw = 80.33 GeV; the systematic uncertainty is 
~ 4 MeV. 



11 



The effect of the resolution is studied by varying the resolution parameter in the pseudo-data, fitting 
to templates with fixed resolution, and collecting the corresponding values of m^,'. This provides the 
relation between the resolution bias and the resulting bias on mw'. 

4^ =0.8 MeV/% 

as illustrated in Figure [5] Injecting the expected precision on the resolution, using the same method 
as above, yields dmw{0() ~ 1 MeV. 

The analysis presented here was originally done in terms of E (rather than transverse energy, Ej) to 
ease comparison with the scale and linearity measurements performed on ATLAS testbeam data [27] . 
In the context of collision data, the analysis can instead be performed in terms of Ej; the energy- 
dependent scale is reconstructed with the same precision as above, in the range 20 < < 70 GeV. 
The propagation to dm]v{ae) and dmw{Oi) is unchanged. 

In addition to the transverse energy dependence, the detector response is in general also a function of 
the lepton pseudorapidity rjf , azimuth 0^ , and time. The physical <pi distributions are however uniform 
in general, and certainly for W and Z events. Hence it is safe to average over <pf, any azimuthal depen- 
dence of the detector scale or resolution then acts as a contribution to the averaged detector resolution. 
Any possible time dependence of the energy response can be treated in the same way, provided the 
analysed W and Z event samples are taken from identical data taking periods ("runs"). It is however 
beneficial to limit the impact of this time dependence on the detector resolution by precisely monitor- 
ing its response as a function of time. 

Although not strictly identical, the rjf distributions in W and Z events are also expected to be very 
similar within the detector acceptance (the difference is below 5% within \rii\ < 2.5, cf. FigurelT]). As 
a first approximation, the same procedure can be applied; the averaging over Tj then assumes that lep- 
tons from W and Z are reconstructed with similar performance, with the same averaging contribution 
to the global detector resolution. The averaging can be improved by reweighting the rjc distribution 
observed in Z events, where the scale is measured, to reproduce the distribution observed in W events 
where the scale applied. The detector response to leptons of given transverse momentum is then iden- 
tical by construction in W and Z decays, up to the statistical precision of the reweighting. As will 
be seen in Section [6l the absolute scale determination is very stable against variations of the under- 
lying physics hypotheses. In particular, it is negligibly affected by PDF uncertainties, which are the 
main factor determining the physical rapidity distribution of the Z boson and its decay products. The 
reweighting does thus not introduce hidden physics uncertainties, and does not affect the discussion 
of other systematic uncertainties. 

The above analysis is performed on the example of the electron channels. As discussed in Section l2!2l 
the muon channels provide similar statistics, and are reconstructed with similar resolution. The present 
results thus equally hold in the electron and muon channels. 

We end this section by noting that other well-known physics probes of detector scale exist, such 
as the low-mass vector resonances and T. An over-constrained scale measurement can also be 
performed by first measuring the ID scale, exploiting muon final states; transporting the ID scale to the 
EMC, using the E/p distribution with isolated electrons; and finally verifying that this indirect EMC 
scale allows to reconstruct unbiased mass peaks for the known resonances in electron final states. 
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Figure 7: Electron (left) and muon (right) rj distributions at reconstruction level, for W and Z events. 



This confrontation of measurements is expected to allow to understand the source of any observed 
non-linearities in terms of magnetic field effects, imperfect alignement, excess of passive material 
in the detector, etc. It will thus be possible to confront several independent probes of the detector 
scale; compatibility between these measurements then validates its use for the measurement of miy. 
This discussion is familiar from the Tevatron mw measurements [9, 10]. The present analysis, using 
Z events only, quantifies the precision achievable at the LHC provided all measurements of the scale 
agree. 

3.2 Lepton reconstruction efficiency 

The observed Jacobian distributions in W events also reflect any pj dependence of the lepton recon- 
struction efficiency. Any difference between the simulation used to produce the templates and the data 
will induce a distortion of the spectrum and cause a bias in the mass fit. 

We again take the electron channel as our main example. The ATLAS electron identification largely 
exploits the shapes of their calorimetric showers [15], which have significant energy dependence. 
Hence, any selection based on these will have a py-dependent efficiency which has to be appropri- 
ately simulated in the templates. Unlike the electrons, no strong pj dependence affects the muon 
reconstruction efficiency. 

3.2.1 Electron efficiency measurements 

Electron reconstruction efficiency can be determined from the data with Z events, using e.g. the so- 
called "tag and probe" method [23], which we briefly summarize here. 

Events are selected with one well-identified electron, and an additional high-T^j-, isolated track. The 
invariant mass of these two objects is required to be within 10 GeV from the nominal Z boson mass. 
Assuming that this selects Z events with enough purity, the identification efficiency is then obtained 
by computing the fraction of events where the second object is indeed identified as an electron. The 
efficiency of the isolation criterion is obtained in a similar way. Simulation studies show that the im- 
pact of backgrounds on the estimation of the efficiency is small compared to the statistical uncertainty. 
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Figure 8: Electron reconstruction efficiency, as determined from fully simulated ee events. 

For the present study we use about 200000 fully simulated Z ^ ee events, from which the efficiency 
is evaluated. The result is shown in Figure [H together with an empirical function describing main 
features of the pj dependence. The following form: 

£(P7-) =£o-aexp(-Zjxpr) (7) 

correctly describes the efficiency in the range relevant for the analysis. 

3.2.2 Propagation to mw^: 5mvi/(£) 

The effect of the efficiency uncertainty is estimated as in the previous section. Template distributions 
are produced at generator level, with varying values of m^y, and applying an efficiency factor accord- 
ing to the best fit efficiency function obtained above. 

One hundred independent pseudo-data samples are generated at a fixed mass {m\N = 80.33 MeV). Ef- 
ficiency functions are applied with parameters drawn randomly within their uncertainties, as obtained 
in the previous section. 

For each sample of pseudo-data, a fit is performed to the W mass. The fitted mass values are his- 
togrammed, and the spread of the histogram gives the corresponding systematic uncertainty. With 
the efficiency determined using 2 x 10^ Z boson decays, the efficiency-induced systematic W mass 
uncertainty is found to be 5mw = 33 MeV. Other functional forms than Eq. |7] yield the same result. 
The most sensitive parameter in Equation |7] is b, the slope in the exponential. It is determined to be 
b = 0.068 ± 0.006, corresponding to a precision of 9%; in other words, dmw/di-eib 4 MeV/%. 

To emphasize the importance of this effect, the same pseudo-data samples ai^e compared to templates 
assuming no /jj -dependence in the lepton reconstruction efficiency (i.e. /{pr) = constant). While 
the same spread is observed, the m(^' distribution indicates an average bias of about 450 MeV. This 
bias vanishes, to first order, when using the /^r-dependent efficiency in the templates. 
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Pt cut < > < 8mi^ > < m^, > < Sni^, > 







(Sref = 1) 


(eref = 1) 


(Sref = f{pT)) 


{£ref = f{PT)) 


Pt 


> 20 GeV 


80.78 


0.033 


80.34 


0.033 


Pt 


> 34 GeV 


80.51 


0.019 


80.34 


0.018 


Pt 


> 37 GeV 


80.44 


0.013 


80.33 


0.012 



Table 2: Average value of and its spread 5m^ , for several lower cuts on the pj range used in the 
mass fit. Numbers are given as obtained from templates assuming a flat efficiency (second and third 
column), and using the efficiency measured in Z events (fourth and fifth column), ni'^'^ = 80.33 GeV. 

Extrapolating to 10 fb^^ i.e. assuming 3 x 10^ measured Z boson decays, an improvement of a 
factor 4 is expected in the efficiency determination. Con^espondingly, we obtain 8m'w{e) ~ 8 MeV. 

3.2.3 Discussion and improvements 

As can be seen in Figure [H the electron efficiency varies most rapidly when pj ~ 20 GeV, and is much 
flatter around the Jacobian edge. Until now, the full pj spectrum, selected as described in Section |2l 
has been used in the mass fits. 

The effect of restricting the lepton pj range used in the fit to higher values is displayed in Table |2l 
Considering the part of the spectrum verifying pj > 34 GeV, for example, reduces 5mw from 33 
MeV to 18 MeV. While avoiding the region with strongest -dependence of the efficiency, the Jaco- 
bian edge is still fully exploited, and the statistical sensitivity is almost unaffected. Extrapolating to 
10 fb^', we obtain a remaining uncertainty of 5m-w{£) ^ 4.5 MeV. 

Note that the results presented here reflect the state of the ATLAS reconstruction software at the time 
of writing this paper. Significantly improved algorithms are described in Ref. [21], notably resulting 
in a smaller /jj -dependence of the electron reconstruction efficiency. The related systematic uncer- 
tainty on my/ should decrease accordingly. The numbers presented here may thus be considered as 
conservative. 

For muons with sufficient momentum to cross the whole detector (p > 6 GeV), no source of ineffi- 
ciency has a strong pr dependence. Hence, the con^esponding induced uncertainty on mw is smaller. 
The above estimate is thus conservative when applied to the muon channel. 

3.3 Recoil scale and resolution 

When using the distribution in the mass fit, pj enters the definition of the observable. This 
quantity, measured experimentally as the vector sum of the transverse energy of all reconstructed 
detector signals (high-pr leptons and low-pj hadronic activity), needs to be precisely described by 
the simulation for the same reasons as above. 

3.3.1 Sensitivity to the recoil scale with Z events 

The W and Z bosons are produced through very similar partonic processes, and thus one expects the 
spectator part of the event (the underlying event) to behave similarly, up to the small phase space 
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Figure 9: Left: statistical sensitivity to as a function of the accumulated Z statistics. Right: 

statistical sensitivity to Op miss . 



difference (mw ^ mz). 

Assuming that the absolute lepton scale and resolution have been measured beforehand (cf. Sec- 
tion l3.1l ). one can measure the recoil scale (a^c) and resolution (Orec) in fully reconstructed Z events, 
where no significant Ef '''^ is expected, by comparing the measured hadronic energy Erec, recoiling 
against the Z boson, to the reconstructed di-lepton four-momentum, pj. Specifically, a^ec and Orec 
are extracted from the peak position and spread of the distribution of Erec/Pj - The results can then 
be used to correct the observed recoil, and hence E'j"^\ in W events. 

Figure |9] shows the expected sensitivity to arec and Orec- With 10 fb^', these pai^ameters can be deter- 
mined with a statistical precision of Sa^c = 5 x 10^^ and SOrec = 6 x 10^"*. 



3.3.2 Propagation to mil/: Sm-w^a^miss), 5myi/{o^miss) 

The effect on m-w is evaluated by systematically varying the recoil scale, producing corresponding 
pseudo-data samples as in the previous sections, and fitting each sample to perfectly calibrated tem- 
plates. We obtain the relation between the bias and the recoil scale and resolution in the form of 
a derivative: 



as illustrated in Figure [TOl Injecting Salmis-: = 5 x 10 ^, we obtain a systematic uncertainty of 
5m-\/v{agmhs) = 1 MeV. Similarly, we find the contribution from the resolution to be 5 
1.5 MeV. These numbers assume that the Z-based calibration can be transported to the W sample 
without additional uncertainty; this is discussed further below. 



3.3.3 Further discussion 

The E'™^'^ calibration can be studied in more detail, using real Z events where one reconstructed lepton 
is artificially removed. In the case of electrons, the removed calorimetric energy should be properly 
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Figure 10: Left: Bias on mw, dm^ = m^' — m'^'" , as a function of the bias on the recoil scale, 5a^,mss. 
Right: 8mw as function of the resolution bias, So^miss. A linear dependence is observed in each case, 
with dm-sff jdreia^mhs = —200 MeV/% and dm-w /dreiO^„ms = —25 MeV/%. 




Figure 11: Left: electron cluster in n Z ^ ee event. Right: the same calorimeter region, after the 
cluster has been removed. The energy in each cell belonging to the electron cluster is replaced by a 
number drawn from a Gaussian with mean and RMS corresponding to detector noise. 



replaced by the expected noise. For muons, also, the minimum-ionizing energy depositions in the 
calorimeters need to be removed and replaced by the expected noise as above. The resulting events 
mimic W events and have a precisely known missing energy, corresponding to the energy of the re- 
moved lepton, which can be compared to the result of the Ej'^^ reconstruction algorithm. 

The lepton removal requires that one can identify and remove the electi^on signal from the struck 
calorimeter cells, while leaving a realistic contribution from noise and hadronic background (see Fig- 
ure [TTll. Several approaches can be tried, such as replacing the contents of the electron cluster cells by 
energy measured away of any high-pj- object in the event (e.g. at 90° in azimuth), or by the average 
expected electronic and hadronic noise. 

To determine the Ej"^^ resolution and possibly correct for biases in its measurement, we consider the 
reconstructed E'™^'^ of Z ^ U events before and after the removal of one lepton, and compare the 
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Figure 12: Transverse view of a) IV ^ ly and \i)Z^ II events. The combined transverse momentum 
of the recoil u, which should match that of the boson, is used to estimate the momentum of the 
undetected neutrino in the W — > decay. The Z boson line of flight is represented, which defines the 
(II , _L) coordinate system. The size of the dotted ellipses represent the resolution on the reconstructed 
objects. 



difference to the transverse momentum of the removed lepton. A non-zero average value of this dif- 
ference points to a bias in the SJ?"* reconstruction. 

Rather than projecting this difference on conventional X and Y axes in the transverse plane, it is best 
to consider the natural frame of the event, with axes parallel (||) and perpendicular (_L) to the Z bo- 
son transverse momentum. Imperfect calibration of the f"™"* reconstruction will show up as biases 
in these distributions, which can then subsequently be corrected for within statistics. The axes are 
illustrated in Figure [121 

This method is tried on a fully simulated sample of Z — > events, with results illustrated in Fig- 
ure [13] As can be seen in this example, a bias is observed in the f"^'" reconstruction along the Z 
line of flight. No bias is observed along the other axes. In this example, the calibration is thus cor- 
rect on average, but the Ej^^^ reconstruction does not respond perfectly to the event-by-event topology. 

As this discussion illustrates, f"^'" reconstruction is a very difficult experimental algorithm to control, 
especially to the level of precision desired here. Therefore, we cannot claim at present that the sensitiv- 
ity quoted in the previous section will indeed be reached. Instead, lacking the proof that the statistical 
enhancement can be fully exploited, we assume an overall uncertainty of 5mvi/(£'™"*) = 5 MeV. This 
number is a factor 3 higher than the purely statistical sensitivity, and a factor three smaller than the sys- 
tematic uncertainty obtained in the recent CDF measurement [10] based on an integrated luminosity 
of 200 pb^' and about 8000 Z events for calibration of the hadronic recoil. 
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Figure 13: Top: resolution on Ef^^'\ projected onto the (X,Y) coordinate system, for unmodified, fully 
simulated Z ^ ee events. Bottom: Ej'^^'^ resolution in the (||,-L) coordinate system. The absence of 
bias along the X and Y axes show that the overall calibration is correct on average, but the observed 
bias along the ||-axis, corresponding to the Z line of flight, indicates imperfect calibration of the 
response to the event-by-event topology. 
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4 Theoretical uncertainties 



We discuss below the uncertainties related to imperfect physics modeling of W production. The 
correlation of the mass measurement with the W width, the impact of final state radiation, and biases 
in the and nij distributions induced by pj and distortions are discussed in turn. 



4.1 boson width: 5mw(riy) 

A change in the W width Fvi/ affects the Jacobian peak, and can cause a bias in the W mass measure- 
ment. To assert the size of this effect, samples with the same W mass but W widths varying in the 
range 1.7 — 2.5 GeVwere produced and subsequently fitted. The relation between Fi^ and mw in the 
fit is linear, with a slope depending on the distribution used in the mass fit. If the W transverse mass 
is used, we find: 

If the lepton transverse momentum is used, we find: 

dmw 



3.2 MeV/% 



1.2MeV/% 



The intrinsic width of the W resonance Fiy has been measured to be 2.141 ±0.041 GeV, while the 
SM prediction is 2.0910 ± 0.0015 GeV [28]. It should be taken into account that the LHC data can 
be expected to improve the precision on the W width as well as on m^. Eaiiier measurements of 
FvK [29, 30] are affected by the same systematic uncertainties as those discussed in this paper. Hence, 
anticipating on our results, we assume that an improvement by a factor five should be achievable, 
respectively leaving 5miy(Fvi') =1.3 and 0.5 MeV for the and pj fits. 



4.2 QED final state radiation: dmw (QED) 

Final state radiation causes significant distortions of the naive, lowest order pj spectrum of the W 
decay leptons. We estimate the stability of the theoretical calculation below, using the PHOTOS pro- 
gram [17] as a benchmai^k. 

The numerical importance of final state radiation is illustrated in Figure [T4j which displays the dis- 
tribution of the measured lepton energy fraction (relative to their energy in the absence of FSR). For 
electrons, measured via calorimetric energy clusters, most of the (coUinearly radiated) photon energy 
is collected in the cluster. The momentum of muons tracks, on the contrary, is measured indepen- 
dently of any photon radiation. The average values of the distributions lie at about 99% of the original 
value, meaning that ignoring the effect entirely would cause a bias on the W mass of about 800 MeV. 
The theoretical stability of the calculation is thus of critical importance. 

In recent versions of PHOTOS, it is possible to switch between several theoretical assumptions. In par- 
ticular, W and Z boson decays can be simulated with photon emission up to 0(a), O(a^), O(a^), or 
with photon emission exponentiation [31]. To study the model differences, we have generated about 
10^ events for each setting, and for each production and decay channel (W — > £v, Z — > M, for i = e,lJ.). 

The average values of the energy fractions discussed above are shown in Figure \T5\ for successive 
theoretical refinements. The different average values for electrons and muons reflect the different 
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Figure 14: Distribution of the measured lepton energy fraction (i.e. relative to their energy in absence 
of FSR). PHOTOS is run in exponentiated mode. The energy of electrons is measured dressed with 
all photon energy radiated within a cone of radius 0.1, corresponding to the size of reconstructed EM 
clusters. Muon momentum is measured bare, after FSR. 



ways their energy or momentum is measured. The calculation appears stable to about l-2x lO^'*, the 
residual differences being compatible with coming from the finite sample statistics only. It is unfor- 
tunately not practical to further increase the samples sizes and quantify the stability to better precision. 

To improve on the above argument, consider the Z boson mass measurement at LEPl [3]. Similarly to 
our case, QED corrections, in the form of initial state radiation off the electron beams, have a large im- 
pact on the Z lineshape, inducing a decrease of the cross-section of about 30%, and a shift of the peak 
position of about 100 MeV. Nevertheless, the theoretical uncertainty on these effects are estimated to 
0.3 MeV, compared to a total measurement uncertainty of 2.1 MeV. The theory of QED radiation thus 
carries negligible uncertainty. 

For the QED induced niiy uncertainty to be as small, the event generators used to produce our tem- 
plates thus need to have similar theoretical accuracy, with the additional complication that the present 
analysis requires an exclusive description of the final state (i.e, a complete description of the photon 
distributions), whereas the Z lineshape analysis only relies on the effective energy of the beams after 
radiation. In Ref. [32], the accuracy of the PHOTOS algorithm is upgraded to NLO accuracy. Simi- 
larly, the HORACE event generator [25] contains QED and weak connections to NLO accuracy. Both 
programs implement photon emission exponentiation. 

We thus assume that ultimately 8mw{QED) < 1 MeV can be reached. This assumption is conditioned 
by the availability of the necessary tools in time for the measurement. 
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Figure 15: Averages of the distributions of Figure[T4l for various PHOTOS settings (see text). 



Let us finally note that W and Z events behave differently under QED radiation. The average energy 
fraction in Z events is 5-7 x 10^^ smaller than in W events, depending on the final state. The energy 
scale measurement (cf. Section 13.11 ) and the W mass measurement should properly account for the 
difference in the respective QED radiation patterns. We will come back to this point in Section [6l 



4.3 W distributions 

The W rapidity and transverse momentum distributions result from the interplay of the proton struc- 
ture functions, and strong interaction effects at the W production vertex. To simphfy the discussion, 
we will consider the longitudinal and transverse distributions independently, as respective results of 
parton distributions and QCD higher orders. 



4.3.1 Rapidity distribution: 5mw{y^) 

The W rapidity distribution is essentially driven by the proton parton density functions (PDFs). Our 
study is based on the CTEQ6.1 structure functions sets [33], which provide, in addition to the global 
best fit, PDFs corresponding to the variation of each diagonal parameter (i.e, the linear combina- 
tion of input parameters that diagonalize the covariance matrix) within its estimated uncertainty. The 
PDF-induced uncertainty for an observable is obtained by computing its value with all sets, taking the 
central value as given by the best fit, and quadratically summing the biases (w.r.t the best fit value) 
obtained from the uncertainty sets. 

As illustrated in Figure [16] (see also Ref. [14]), the current PDF uncertainties induce an uncertainty 
in the W rapidity distributions which, thi^ough acceptance effects, propagates a systematic uncertainty 
on the W mass determination of ~25 MeV. We present below an attempt to estimate how this will 



22 




Figure 16: Bias on mw obtained when varying the proton PDFs within their uncertainties. Each point 
on the abscissa coneponds to a given PDF set: set is the best fit, and gives bias up to the statistical 
uncertainty of the fit; sets 1-40 are the uncertainty sets, each inducing a given bias on niw- The total 
uncertainty on mw is given by the quadratic sum of the biases, giving 5mw ~ 25 MeV. 



improve with the LHC data. 

At the LHC, W and Z particles are essentially produced through sea quark interactions; the influence 
of valence quarks is small. Low-x, high-2^ sea quarks mainly evolve from higher x, lower glu- 
ons, and a consequence from perturbative QCD flavour symmetry is that up to initial asymmetries 
and heavy-quark mass effects, the different quark flavours should be represented democratically. This 
then implies that the impact of sea quark PDF uncertainties on W and Z production should be very 
similar. In other words, when varying PDFs within their uncertainties, one expects a strong correlation 
between the induced variations of the W and Z distributions. 

This is confirmed by Figure [iTlEl On the left, the coiTclation between the widths of the W and Z 
boson rapidity distributions is displayed. We choose to use the distributions RMS, denoted and 
r^, to quantify their width. The current CTEQ6. 1 prediction, = 2. 16 ±0.03, will be refined to a 
precision of 5r^ = 0.001. Exploiting Figure [T7] (right), which quantifies the correlation between 
and r^, this can be translated into a prediction of the W boson rapidity distribution, 5r^ = 0.0013, to 
be compared to the current prediction = 2.24 it 0.03. 

One thus expects an improvement on the Z rapidity distribution by a factor ^^30. This is also illustrated 
in Figure [181 where two extreme predictions (with cunent knowledge) of the Z rapidity distribution 
are compared with an example distribution representing the same measurement. Given the residual 
decorrelation between the W and Z distributions, this translates into an improvement on the W rapidity 
distribution by a factor ~ 23. 

Starting with dmw{yw) ~25 MeV, putting in a precise measurement of the Z rapidity distribution at 

'^This plot is reminiscent of Figure 2 in [34], displaying similar correlations in the production rates. Note that for our 
purpose, normalizations are irrelevant and we are interested only in the distributions. 
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Figure 17: Left : conxlation between the spreads (RMS) r|f and of the W and Z rapidity distribu- 
tions, when varying the CTEQ6.1 PDFs within their estimated uncertainties. The fitted pseudo-data 
are scaled to an integrated luminosity of 10 fb^ ^ Right : distribution of the ratio /ff, again varying 
the PDFs within their uncertainties. The spread of the ratio distribution is 4 x 10^^. 




Figure 18: The line histograms represent two extreme predictions for the Z rapidity distribution, as 
given by the CTEQ6. 1 PDF sets. The points are pseudo-data, obtained with the central set, and scaled 
to an integrated luminosity of 10 fb^^. 
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Figure 19: Upper plot: the outer histogram represents the complete rapidity distribution for W produc- 
tion at the LHC; the inner histogram represents the range selected by the condition \ri(\ < 2.5. Lower 
plot: the outer histogram represents the complete rapidity distribution for Z events. The innermost 
histogram is obtained requiring two decay leptons within |t7£| < 2.5; the intermediate histogram is ob- 
tained when allowing one electron within |t7/ | < 4.9. The two symmetric histograms at high rapidity 
correspond to the LHCb muon acceptance. 



the LHC, and exploiting the strong correlation between the W and Z production mechanisms, we thus 
anticipate a final uncertainty from the description of the W rapidity distribution of 5mw{yw) ~1 MeV. 

In practice, the analysis will of course proceed via a formal QCD analysis to the LHC data: the mea- 
sured Z differential cross-section da/dy, together with other measurements (see below), will be fed to 
parton distribution fits, and the systematic 5my\r{y^) from the improved PDF sets will be evaluated as 
above. The present discussion however allows to estimate the expected improvement while avoiding 
these complications. 

Let us also note that Z rapidity distribution can be analyzed over a domain that fully includes the range 
relevant for W production. In ATLAS (as in CMS), the usual Z acceptance, given by \ri(\ < 2.5 for 
both decay leptons, can be extended in the electron channel by allowing one of the electrons to be de- 
tected within |T]e| <~ 4.9. In addition, high-rapidity Z events will be produced and detected at LHCb 
(for example, the geometric acceptance of the muon detector is approximately 2.1 < |t]^| < 4.8). 
Accounting for this, and as illustrated in Figure [l9l the W rapidity range selected for the niw measure- 
ment is entirely included in the Z one. This remains true in terms of the parton momentum fractions. 

We conclude this section with some caveats. The above results partly are a consequence of the as- 
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sumed flavour and charge symmetry in the low-;c proton; notably, the parton parametrisations used in 
the fits used above assume that d{x) = d{x) = u{x) = u{x) at low-;c, and 5 = 5 at all x. This implies 
the strong correlation discussed above, since the Z production rate is proportional to uu + dd+ . . ., 

and the W rate is proportional to ud + dU + It is thus important to quantify the dependence of our 

result on these hypotheses. 

The anti-quark flavour asymmetry m — </ was measured to be non-0 in the region 0.015<;c<0.35, and 
^ 50 GeV^ [35, 36], in contradiction with the flavour symmetry assumption. The relative asym- 
metry, {a — d) / {{a + d) , is however of the order ~ 10^^, decreasing towai^ds higher Q^. Starting from 
u = d and full correlation between W and Z production (i.e. W and Z distributions have the same 
rate of change under PDF variations), d induces a decorrelation of order {u — d)/ {11 + d) x [u — 
d) j (u+d), where both factors are of order 10^^ (see for example Figure 1 in [33]). Hence, even in the 
presence of non- vanishing u — d, the freedom of the W distributions is very limited once Z ones have 
been precisely measured. We thus assume that our estimates remain correct; nevertheless, measure- 
ments of the W charge asymmetry, sensitive to ii — d, will allow to verify this hypothesis. Additional 
information will be provided by measuring mw in W+ and events separately. 

The proton strangeness asymmetry, s{x) —s{x), is constrained by neutrino scattering data [37-39]. 
The relative asymmetry is rather small, even at low Q^: {s — s)/{s + s) ~ 10^^ at = 10 GeV^. It 
will only become smaller at ~ ni^, where most of the strange sea is generated radiatively. We 
consider, as above, that the contribution of the asymmetry is small in terms of the overall W pro- 
duction and its uncertainty. However, the impact on the mw measurement would need to be studied 
specifically. At the LHC, the analysis of + c/c production should provide additional insight. 

Finally, one may argue that the influence of heavy quark PDFs on W and Z production is different, thus 
a source of decorrelation between the two processes. The charm quark contribution to W production 
is significant (~ {VcsCS + Vcdcd + c.c.)), but smaller for Z production (~ cc). On the other hand, the 
b-quaik content contributes to Z production (~ bb), but negligibly to W production (~ {Vctcb + c.c.)), 
due to the smallness of the off-diagonal third generation CKM matrix elements. These differences are 
however accounted for by the present analysis, since the heavy quark PDFs aie included the CTEQ6.1 
PDF sets; heavy flavours are actually understood to cause in part the small decorrelation between the 
W and Z boson distributions. Our conclusions thus remain unchanged. 

The present study has been repeated using the MRST2001 PDF sets [40]. The same con^elation is 
observed between r|f and r^, and the same result is obtained. Non-global parton density fits, such 
as those performed by the HI and Zeus experiments, are based on similar hypotheses and claim 
slightly smaller uncertainties [41], again preserving our result. Finally, during the course of this work, 
CTEQ6.5 PDF sets became available [42], which improves on the treatment of heavy quai^k masses 
in the QCD evolution. The flavour symmetry assumptions are however unchanged, so that the present 
discussion is not affected. 

4.3.2 Transverse momentum distribution: SniwiPx) 

The prediction of vector boson pj distributions at hadron colliders has long been an active sub- 
ject [24,43,44]. It is also a crucial input for the W mass analysis, especially when using the pj 
observable. We discuss below the impact of p^ uncertainties on the W mass determination in this 
hypothesis. 
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The measurable pj and /^f distributions are the resuh of several effects, most notably the repeated, 
partly non-perturbative parton radiation occurring in the transition from the low-Q^ proton towards the 
hard process (commonly refeiTcd to as parton showers, or soft gluon resummation). Another source is 
the transverse momentum intrinsic to the partons in the proton. We choose not to discuss these effects 
separately. Rather, reckoning that although W and Z production differ in several respects (the coupling 
to initial partons is different in both phase space and flavour), the non-perturbative mechanisms are 
universal, we evaluate how precisely their combined effect can be measured in neutral cun^ent events, 
and how this improves the IV predictions. Notice that heavy flavour PDF have caused only a small 
decorrelation between and Z events in the previous section; this is assumed to remain true in this 
discussion. 

First, the relation between the bias in the modeling of and the measurement of my^ is investigated 
by applying scaling factors to the distributions in our pseudo-data, deducing the corresponding 
Pj distributions, and fitting mw against un-distorted templates. The bias in m^/^ appears to be a lin- 
ear function of the pj mis-modelling, with a slope of order 0.3, meaning a 3 MeV bias on p^ results 
in a 1 MeV bias on miv, when exploiting the pj distribution. When is used, the effect is negligible. 

Neutral current dilepton events allow to measure the pj distribution, as a function of mass, over 
a large mass range. Assuming usual selections, this distribution will be measured precisely for 
30 < Mu <~ 200 GeV. This large lever arm, in addition to the very precise determination of the 
Pj distribution on the Z peak, provides a precise control of do jdpj when M(i ^ rtiw. This is illus- 
trated in Figure |20l which displays the dilepton mass dependence of its average transverse momentum, 
<Pj >, as predicted by PYTHIA. 

On the Z peak, pj will be known to about 7 MeV with an integrated luminosity of 10 fb^'. Thanks 
to the Drell-Yan continuum, the accuracy in the region of mw is still ~8 MeV. This precision can be 
used to constrain the non-perturbative parameters governing the parton shower or resummation com- 
putations, and to predict the p'j distribution with similar accuracy. This leads to an uncertainty on 
of about 3 MeV. 

Arguably, the p^ distribution cannot be summarized by its mean value. However, in the low p^ 
region (selected by the recoil cut, cf. Section |2ll, it can be empirically described by a two-parameter 
function. As an exercise, the mass-dependence of the pai^ameters were determined on Drell-Yan 
events, their values and uncertainties in the miy region were used to produce pj pseudo-data as above, 
and corresponding fits to m^/v were performed. The spread in myv resulting from the uncertainty in the 
empirical parameters was found compatible with the above estimate. 

5 Environmental uncertainties 
5.1 Backgrounds 

The leptonic W final states benefit from low backgrounds, mostly coming from vector boson de- 
cays; notably W — > t(— > £vv)v (in^educible), Z ^ II (where one lepton is not reconstructed), and 
Z — > t(— > ^vv)t. QCD dijet events will, despite their large cross section, not be dominant. The 
backgrounds from tt and W+W^ events are negligible. The systematic error on m-^ arises from un- 
certainties on the background shape and normalization in the fitting range of the and spectra. 
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Figure 20: Top : Dilepton invariant mass spectrum, from inclusive neutral current events (7 and Z 
exchange are included). Bottom : dilepton average pT as a function of the dilepton invariant mass. 
The W-mass region is strongly constrained by the lever arm provided by the Z peak and the Drell- 
Yan rise at low mass (note the improved precision in these regions). The points correspond to a 
measurement with 10 fb^'. 
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Figure 21: Signal and backgrounds in the pj distributions, for W ^ ev (left) and W jiv (right). 
The histograms correspond, from bottom to top, to Z — > tt, Z — > W — > tv and W — > iv. 



Uncertainties on the W and Z background size, relative to the signal size, depend on cross-sections, 
branching fractions and acceptances. These aie obtained from the PDG [28] and take into account the 
studies described in Section \i?2\ and 1431 Note that in contrast to the studies presented until now, the 
background uncertainty does not scale with statistics. 

The background shapes are determined from simulation. They aie essentially unaffected by variations 
in the production, decay, and resolution model, and play only a minor role in the overall systematic 
errors. For QCD background, as a separate study, both normalization and shape will have to be mea- 
sured directly from the data. The pj distributions, including signal and backgrounds, are illustrated 
in Figure |2T] 



W — > TV events: The largest background is from W — > TV events, where the t decays into a lepton. 
This background is iiTcducible, as the final state is identical to the signal; however, its pj and ni^ are 
on average lower, leaving a tail into the fitting range. Though being the main background, its uncer- 
tainty is small, as only T decay parameters and the acceptance enter, with respective uncertainties of 
1% and 2.5%. 



Z U events: The second largest background is from Z ^ II events, where one lepton is either 
undetected or not identified. This background can be reduced using a Z veto rejecting events, where 
the lepton and a second isolated object (track and/or cluster) form an object with an invariant mass 
between 80 and 100 GeV (see Figure l22l) . Due to the high mass of the Z boson, the pj distribution 
extends well into the fitting range. The distribution is again at low values, due to the smallness 
of missing momentum. The size of this background has uncertainties from both the W to Z cross 
section ratio Rwz, and from the acceptance/veto efficiency. It is expected to be larger for muons than 
for electrons, as the former cannot be vetoed for |t] | > 2.7. 
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Invariant ee mass (GeV) 

Figure 22: Distribution of invariant mass between lepton and a second isolated object (track and/or 
cluster) in Z — > events where only one lepton is identified. Events in the range 80-100 GeV are 
rejected. 



Z ^ TT events: A small background originates from the Z — > TT process, where one T decays lep- 
tonically, while the other is not identified. While the cross section for such a process is small, it 
contains significant fi™'". 



Jet production: The QCD background cannot be obtained reliably from simulation. It will thus 
have to be measured directly from data. For the Run I W mass measurement at CDF, this background 
could be estimated to a precision of ~50% [45], limited by lepton identification performances and 
statistics. At ATLAS, a precision of ~10% is assumed in the electron channel, where this background 
is expected to be significant. The assumed improvement is justified by the superior granularity and 
resolution of the EM calorimeter [15]. The muon final state is less contaminated by jet events, muons 
being measured behind all calorimetry. A specific background is however constituted by muons from 
hadron decays in flight. As we have no measure of the uncertainty on this background, our results 
implicitly assume it is small. We stress that these estimates are essentially quantitative. A realistic 
estimate of their impact on the measurement will only be possible with data. 



Overall impact: We now estimate the overall impact of the backgrounds. The background shapes 
can be empirically described by an exponential function in the fitting range, as illustrated in Figure [23] 
on the example of the W ^ TV background. The systematic uncertainty on is then derived by 
varying the function parameters within their uncertainties as estimated above. The systematics uncer- 
tainty induced by the background shapes amounts to 20% of that induced by the normalizations. 

The overall effect is obtained by repeating this procedure for all backgrounds. Table [3] summarizes 
background uncertainty and its impact on the W mass determination. 
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dis1_plej 



dist_pte 


Entries 


10000 


Mean 


30.21 


RMS 


4.243 


1 ndf 


3575 / 37 


pO 


2.397e+04 + 8383 


p1 


18.5 + 2.4 


p2 


8.565 ± 0.583 




Figure 23: W — > TV background shape in the fitting range (indicated by dashed Unes) for irvf (top) 
and pj (bottom). 



Background 


Variable 


Error 


Derivative 


Impact (MeV) 


W ^zv 


W 
rrij 


2.5 % 


-0.5 MeV/% 


1.5 




Pt 


2.5 % 


-0.7 MeV/% 


2.0 




W 
nij 


2.8 % 


0.08 MeV/% 


0.22 




Pt 


2.8 % 


0.09 MeV/% 


0.26 




W 


4.5 % 


0.02 MeV/% 


0.09 




Pt 


4.5 % 


0.03 MeV/% 


0.14 


QCD events 


w 
nij 


10 % 


0.04 MeV/% 


0.40 




Pt 


10 % 


0.05 MeV/% 


0.50 


Total 


W 
Pt 






1.6 
2.1 



Table 3: Table of backgrounds along with its uncertainty, derivative, and impact on niw. The overall 
systematic uncertainty from backgrounds is about 2 MeV. 
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Combining the systematic errors from the backgrounds yields a total of 1.6 and 2.1 MeV for the iri^ 
and pj distributions, respectively. 

5.2 Pileup and underlying event 

The soft hadronic activity accompanying the hard process (underlying event), and the overlap with 
soft events produced in the same bunch crossing (pile-up) generate additional particles that contribute 
to the detector occupancy. In particular, the additional calorimetric energy overlaps with the electron 
signal and distorts the electron scale measurement. 

Typically, a soft event produces about 10 particles per unit rapidity (integrated over (^), with average 
transverse momentum pj ~ 500 MeV [46,47]. An electron cluster of typical size 5t] x 5^ ~ 0.1 x 0.1 
is expected to contain about 40 MeV of hadronic background, to be subtracted from the electron signal. 

In particular, the hadronic background may have a non-negligible 2^-dependence, generating a non- 
universality between W and Z events. These effects are small but need to be properly accounted for 
when aiming at a precision on the absolute electron scale of5a/a~2x 10^^. 

This aspect was not studied here, but we follow the argument of [15]. By measuring the energy flow 
away from any high-7?robjects, as a function of T], independently in W and Z events, a 2% precision 
on the hadronic energy flow looks achievable. Such a result would bring down the size of the effect 
from 40 MeV to about 1 MeV. 

We thus conclude that although soft hadronic interactions generate shifts in the energy measurements 
that are large compared to the statistical sensitivity to mvi/, these shifts can be measured in the data 
with sufficient accuracy. The final contribution to Sm^y is small. 

This source of uncertainty affects the electron scale; the muon scale is not affected. The impact on the 
recoil measurement is not discussed here; this section is thus relevant for pj based measurements. 

5.3 Beam crossing angle 

At the LHC, the proton beams are brought to collision at a crossing angle of 142.5jUrad [48]. In terms 
of momentum, this translates into a 7000 GeV x 142.5 x 10^^ 1 GeV boost in the horizontal plane 
(.x-direction), per beam proton. However, in the simulation protons collide head-on, giving rise to a 
systematic shift in p^ of all particles produced. 

Figure |24] shows the difference in the transverse W momentum before and after taking this effect into 
account, b.p^ = - Px^""', which is expected to be up to Ap^ = mw ■ 142.5 x 10"^ « 11 MeV. 
However, since the W boson line of flight has azimuthal symmetry, the impact on the W transverse 
momentum distribution is smaller, as most of the effect is averaged out by the rotational symmetry. 

The size of the effect is estimated as usual, by including the p^""'" in the pseudo-data and letting the 
templates unchanged. We find that the effect is smaller than 0. 1 MeV. 
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Figure 24: Distribution of difference in the transverse W momentum resulting from the boost, Ap 
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6 Correlations 

So far, all main sources of systematic uncertainties have been investigated independently. Before we 
move to the combination of our results, we need to address the question whether important congela- 
tions are to be expected between the sources. It is, however, beyond the scope of this work to discuss 
this issue extensively, and we limit this section to the most important examples. 



The uncertainty related to the absolute scale has the strongest lever arm on the determination of ntw 
(5mw/5a = 1). Therefore, we investigate below whether uncertainties which affect the W mass 
measurement can also bias the absolute scale. 



6.1 Absolute scale vs. lepton reconstruction efficiency 

We repeat the procedure described in Section 13.11 As before, Z boson invariant mass templates are 
produced for different scale and resolution hypotheses, and pseudo-data with scale parameters to be 
determined. The impact of a -dependent lepton reconstruction efficiency is assessed by assuming 
perfect efficiency in the templates (e = 1), and injecting the efficiency function discussed in Sec- 
tion [3]2]in the pseudo-data. 

The result is displayed in Figure |25] It appears that the injected inefficiencies merely induce a re- 
duction of statistics, and hence some loss of precision in the scale determination, but no appreciable 
bias: in spite of the reduction in statistics, the reference invariant mass distribution is not signifi- 
cantly distorted. Note that, since the efficiency is assumed perfect in the templates, and realistic in 
the pseudo-data, any observed bias would have been a large overestimation of the effect, representing 
100% uncertainty on the effect. 



6.2 Absolute scale vs. PDFs 

Similarly as above, and also as in Section 14.3.11 we compare Z boson mass templates produced with 
the CTEQ6.1 central set to pseudo-data produced with the 40 uncertainty sets. 
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Figure 25: Scale determination using Z mass templates assuming perfect identification efficiency. 
Inner pai^abola: perfect efficiency is also assumed in the pseudo-data; outer parabola: the pseudo-data 
incorporate a p 7- -dependent efficiency. 



The results of the 40 fits are displayed in Figure [26l in the form of biases with respect to position of 
the mass peak obtained in the templates. The CTEQ6.1 uncertainty sets induce typical biases of ~0.5 
MeVwith respect to the central value. Summing over all uncertainty sets gives a total scale uncertainty 
of about 2.5 MeV. This translates into 5mw '-^ 2.2 MeV. 

In other words, with current knowledge, the PDF uncertainties induce a direct systematic uncertainty 
of about 25 MeVv/a distortions of the W distributions (cf. Section l4.3.1l) . and an indirect uncertainty 
of 2.2 MeVv/<3 distortions of the Z lineshape, propagating to the absolute scale determination. 

Hence, the conclusions of Section 14.3.11 are essentially unchanged. Using measurements of the Z 
boson distributions, the PDF induced systematic uncertainty should drop to about 1 MeV. 

6.3 Absolute scale vs. QED corrections 

QED corrections affect the determination of the absolute scale in two ways. First, as was mentioned 
in Section l4!2l the observed W and Z decay lepton spectra are strongly affected by photon emission. 
This effect needs to be taken into account properly when producing the Z mass templates. 

In muon final states, the theoretical distributions are based on the final muons, after simulation of 
the QED photon emissions. Final state electrons cannot be separated experimentally from the mostly 
coUinear photons. Hence, the simulation needs to reproduce this recombination precisely. This de- 
mands precise theoretical control of the photon distributions, an aspect which seems under sufficient 
control (cf. Section l4!2l ). Likewise, a precise description of the detector geometry and EMC shower 
development in the simulation are needed to properly simulate the fraction of photon energy recom- 
bined in a given electron cluster. 

Secondly, as a consequence of the above, the absolute scale extracted from Z events actually corre- 
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Figure 26: Bias on mz obtained when varying the proton PDFs within their uncertainties. Each point 
on the abscissa correponds to a given PDF set: set is the best fit, and gives bias by definition; sets 
1-40 are the uncertainty sets, each inducing a given bias on mz- The total uncertainty is given by the 
quadratic sum of the biases, giving 8mz ^ 2.5 MeV. 



sponds to a mixture of photons and electrons. In ATLAS, the EMC response to electrons and photons 
is different by about 1%, an effect coming from calorimeter geometry (because their showers develop 
differently, electrons and photons of a given energy do not "feel" the same sampling fraction) and 
from the passive material in front of the EMC, which causes early showers or conversions, with dif- 
ferent probabilities for both particle types [49]. It is thus important to know whether W and Z behave 
similarly in this respect, and if any difference is well understood theoretically. 

As is shown in Figure |22l the electron energy fraction in EM clusters differs by about 0.6% between 
W and Z events, meaning that the energy scale measured in Z events needs to be corrected by a factor 
1% X 0.6% = 6 X 10^^. Failing to take this factor into account would induce a bias of 5 MeVon the 
mw fit. However, Figure [27] also shows a good stability of the theoretical prediction. Hence, although 
this correction is not negligible, it does not carry a significant uncertainty. 



7 Impact on the W mass measurement 

We summarize below our main results. Table |4] recalls the main systematic contributions to the pj- 
and m^^ -based wiiy measurement, with 10 fb^' of data. In both tables, numbers are given for the elec- 
tron and muon channels separately when applicable. 

The major difficulty is, as expected, the determination of the absolute energy scale of the final state 
leptons and the hadronic recoil. The analysis of the Z peak however allows to strongly constrain the 
lepton scale uncertainty. The analysis is non trivial, because in addition to the Z mass parameters, 
many other effects enter the theoretical description of the lineshape; most notably, QED radiation. 
Although the effect is lai^ge, the theoretical understanding is adequate, as the LEPl Z mass measure- 
ment indicates. The Z mass measurement relies on an analytical formulation of the inclusive radiation 
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R,(W) - R^(Z) 




PHOTOS1 PHOTOS2 PHOTOS3 PHOTOS Exp 




PHOTOS1 PHOTOS2 PHOTOS3 PHOTOS Exp 



Figure 27: For electron final states in W and Z events, the energy fraction deposited by electrons 
in reconstructed electromagnetic clusters (l-Re is photon energy), for various PHOTOS settings (see 
Section 1121). 



spectrum; the W mass measurement at the LHC however requires a complete Monte-Carlo implemen- 
tation, providing an exclusive description of the final state at the same level of precision. Such tools 
are critically needed in the context of this measurement. 

The analysis of the transverse mass requires in addition a precise calibration of the hadronic recoil 
using Z events, and an unbiased transport of the calibration to W events. Such an algorithm is not 
discussed here; the corresponding systematic uncertainty assumed here is a compromise between the 
high statistical sensitivity of the in situ calibration in ATLAS, and the actual result recently obtained 
at the Tevatron [10]. 

The electron channel appears somewhat more difficult than the muon channel. The first reason is the 
pr-dependent electron identification efficiency, which distorts the Jacobian distributions; this effect is 
essentially absent in the muon channel. The second reason is again related to QED radiation: since the 
muons do not recombine with the emitted photons, the description of the effect is purely theoretical. 
In the case of electrons, a large fraction of the radiated energy is included in the electron cluster. De- 
termining this fraction requires a precise description of the detector geometry and reliable simulation 
of EM showers. 

We estimate that uncertainties related to the description of the and pj distributions will be small 
once the Z differential cross-section will have been measured. As discussed in Sections r4.3.1l and[ 4.3.2[ 
this result relies on two assumptions. The first assumption concerns the light quark flavour and charge 
symmetry in the low-x, high-Q^ proton. We estimated that relaxing these hypotheses within bounds 
allowed by the existing data is unlikely to invalidate our result. Another assumption is that the non- 
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perturbative mechanisms controling the and distributions remain essentially universal, although 
heavy flavour contributions to W and Z production are different. The effect of heavy flavours on the 
distribution has been studied [50], but a study comparing these effects on W and Z production is 
cuiTcntly lacking. The largest remaining systematic comes from the modeling of pj , in the pf. -based 
measurement, contributing a 3 MeV uncertainty. The -based measurement is more stable in this 
respect, but suffers additional experimental complications related to the experimental control of the 
^miss reconstruction. 

Backgrounds contribute an uncertainty 5m^)^ ^ 2 MeV. Of all components, the background from jet 
production is the least well known, but its contribution is expected to be small. We did not investigate 
the possible impact of cosmic rays and hadron decays in flight, which occur in the muon channels, but 
Tevatron experience indicates the impact is small. 

All in all, a total uncertainty of about 7 MeV can be achieved, in each channel, using either the pj 
or the nij method, with the equivalent of 10 fb^' of data. Most sources of systematic uncertainty 
seem to scale with the accumulated Z statistics; notable exceptions are backgrounds, QED radiative 
corrections and the underlying event. Their contribution to 8m-w is however subdominant. Combining 
channels, and allowing for more data, we can therefore expect further improvement. 

Let us briefly compare our results with the recent prospects presented by the CMS Collaboration [14]. 
We base our comparison on the p^-based my/ measurement and 10 fb^' of data. CMS claims 2 MeV 
from the absolute scale, agreeing with our average scale result of Section !?. 1.11 A simplified treatment 
of non-Unearities leaves a systematic uncertainty of 10 MeV, and the assumed 8% relative knowledge 
on the resolution contributes 5 MeV; these numbers can be compared to the 4 MeV we obtain in 
Section [3. 1.21 We include a discussion of the reconstruction efficiency uncertainty, which is omitted 
in [14]. On the theoretical side, the present note and Reference [14] agree on the initial uncertainties 
related to PDFs and the description of the W transverse momentum distribution. Our improvements in 
this respect rely on an analysis of the constraints provided by the analysis of the Z boson differential 
cross-section at the LHC. Finally, we claim a statistical sensitivity of about 2 MeV, compared to 
15 MeV in [14]. This is explained by CMS choosing to base the W templates on measured Z events 
{via the scaled observable method, or scaling the kinematics event by event), thus paying for the 
smaller Z boson production rate. Such a procedure is in principle justified by the reduction of other 
systematic uncertainties, but as we saw throughout this paper this does not seem to be a worthy trade. 

8 Conclusions and perspectives 

We investigated the most important systematic uncertainties affecting the W mass determination at 
the LHC, and found that the analysis of Z production constrains the systematic uncertainties to a total 
of about 7 MeV per channel, exploiting 10 fb^' of data. Combining independent measurements may 
bring further improvement. 

Among all investigated sources of systematic uncertainty, two items in particular rely on assumptions. 
The first one concerns the treatment of QED radiation. We argued that the theory is under very good 
control, having notably allowed a very precise Z mass measurement at LEPl, where QED effects aie 
large, but the uncertainties finally have an almost negligible contribution. To preserve this situation 
at the LHC, the m-w measurement requires QED simulation tools providing the same level of accuracy. 
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Source 


Effect 


dm-w/dreia (MeV/%) 


5reia (%) 


5mv/ (MeV) 


Prod. Model 


W width 


1.2 


0.4 


0.5 




distribution 


— 


— 


1 




Pj distribution 


— 


— 


3 




QED radiation 


— 


— 


<1(*) 


Lepton measurement 


Scale & lin. 


800 


0.005 


4 




Resolution 


1 


1.0 


1 




Efficiency 


— 


— 


4.5 (e);<l(At) 


Recoil measurement 


Scale 


— 


— 


— 




Resolution 


— 


— 


— 


Backgrounds 


W ^zv 


0.15 


2.5 


2.0 




Z 


0.08 


2.8 


0.3 




Z ^ TT 


0.03 


4.5 


0.1 




Jet events 


0.05 


10 


0.5 


Pile-up and U.E 








<l(e);~0(At) 


Beam crossing angle 








<0.1 


Total {p'-j) 








-7 (e); 6 (At) 




Source 


Effect 


dmw/dreia (MeV/%) 


5reia (%) 


5mw (MeV) 


Prod. Model 


W width 


3.2 


0.4 


1.3 




distribution 






1 




distribution 






1 




QED radiation 






<1 (*) 


Lepton measurement 


Scale & lin. 


800 


0.005 


4 




Resolution 


1 


1.0 


1 




Efficiency 






4.5 (e);<l(At) 


Recoil measurement 


Scale 


-200 








Resolution 


-25 








Combined 






5 (**) 


Backgrounds 


W^rv 


0.11 


2.5 


1.5 




Z 


-0.01 


2.8 


0.2 




Z^zz 


0.01 


4.5 


0.1 




Jet events 


0.04 


10 


0.4 


Pile-up and U.E 








<l(e);~0(M) 


Beam crossing angle 








<0.1 


Total (m^) 








~8 (e); 7(At) 



Table 4: Breakdown of systematic uncertainties affecting the m-w measurement, when using the pj 
distribution (top) and the distribution (bottom). The projected values of SreiOi ai^e given for a 
single channel and assume an integrated luminosity of 10 fb^^. The QED induced uncertainty (*) is 
realistic given the precision claimed for the Z boson mass measurement at LEPl, but assumes that the 
needed theoretical tools will be implemented in time for the measurement. The recoil measurement 
uncertainty (**) has not explicitly been quantified here, but is conservatively extrapolated from recent 
Tevation experience. See text for discussion. 
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The second assumption concerns the effect of the Ught and heavy flavours in the proton. Releasing 
the light flavour symmetry assumption in use in the current global QCD fits will cause a decorrelation 
between W and Z production at the LHC. This decon^elation can be expected to be small, but will have 
to be measured at the LHC, notably using the rapidity-dependent W charge asymmetry and the study 
of associated W/Z + charm production. Similarly, heavy flavour PDFs generate some decorrelation. 
This decorrelation was verified to be small in the and distributions, and the same was assumed 
true for the and pj distributions. To verify this assumption requires a theoretical study comparing 
the heavy flavours influence on soft gluon resummation in W and Z events. 

A number of sources have not been studied explicitly, notably the recoil measurement, affecting the 
distribution; the underlying event, affecting the electron energy scale; and W polarization effects, 
affecting the leptonic angular distributions. Other sources, like backgrounds from jets, cosmic muons, 
or induced by the machine can only be studied reliably using real data. We believe these mechanisms 
can be brought under sufficient control, on the time scale of the LHC measurement of m^- 

The results presented here have only exploited Z boson measurements. Many other calibration pro- 
cesses exist, that give additional constraints on the detector performance and on the physics mecha- 
nisms influencing W production. While first providing a way to verify the robustness of the Z-based 
calibrations, these processes can help to reduce the uncertainties further in the case of consistent 
results. We reserve these refinements to the analysis of the forthcoming LHC data. 
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